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DISCLAIMER 



Although each program has been tested by its contributor, 
no warranty, express or implied, is made by the contributor 
or any User's Group, as to the accuracy and functioning of 
the program and related program material, nor shall the fact 
of distribution constitute any such warranty, and no responsibility 
is assumed by the contributor or any User's Group, in connection 
therewith. 



COMMON USEES GROUP PROGRAM REVIEW AND EVALUATION 
(fill out in typewriter, ink or pencil) 

Program No. Date 

Program Name: _ 

1. Does the abstract adequately describe what the program is and what Yes No_ 

it does ? 

Comm ent 

2. Does the program do what the abstract says ? Yes No_ 

Commen t 

3. Is the description clear, understandable, and adequate? Yes No_ 

Comment 

4. Are the Operating Instructions understandable and in sufficient detail? Yes No_ 

Comm ent 

Are the Sense Switch options adequately described (if applicable) ? Yes No_ 

Are the mnemonic labels identified or sufficiently understandable ? Ye s No_ 
Comment 



Q 



5. Does the source program compile satisfactorily (if applicable)? Yes No_ 

Comment 



Does the object program run satisfactorily? Yes No_ 

Comment 



7. Number of test cases run . Are any restrictions as to data, 

size, range, etc. covered adequately in description? Yes No_ 

Comment 

8. Does the Program meet the minimal standards of COMMON? Yes No_ 

Comment 

9. Were all necessary parts of the program received? Yes Nq_ 

Comment 

10. Please list on the back any suggestions to improve the usefulness of the program. 
These will be passed onto the author for his consideration. 

Please return to: Your Name 

Company 



Mr. Richard L. Pratt 

Data Corporation Address 
7500 Old Xenia Pike 



Dayton, Ohio 45432 Users Group Code _____ 

O : 

THIS REVIEW FORM IS PART OF THE COMMON ORGANIZATION'S PROGRAM REVIEW AND. 
EVALUATION PROCEDURE. NONME MBERS ARE CORDIALLY INVITED TO PARTICIPATE 
IN THIS EVALUATION. 
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Computation of Equilibrium Composition 
and Temperature of Chemical Reactions 

By 

Mae E. Meads 
and 

George E. McGowan 



Direct Inquiries to: 



Mae E. Meads 

Engineering Computer Applications Group 
The Baltimore Gas and Electric Company 

Baltimore, Maryland 21203 

1620 Users Group No. 123$ 



March 9,. 1966 



Modifications or revisions to this program, as they occur, 
will be announced in the appropriate Catalog of Programs 
for IBM Data Processing Systems. When such an announce- 
ment occurs, users should order a complete new program 
from the Program Information Department. 




1620 USERS GROUP LIBRARY 
PROGRAM ABSTRACT 

1. TITLE (If subroutine, state in Title): Computation of Equilibrium Composition and 

Temperature of Chemical Reactions 

2. Author; Organization: Mae E. Meads, George E. McGowan - Baltimore Gas and Electric 

Company. Subject Classification: 

Date: 1 November. 196$ Users Group Membership Code: 123$ 

3. Direct Inquiries to Name: Engineering Computer Applications Group. Baltimore Gas and 
Electric Company. Baltimore. Maryland 21203 Phone: $39-8000 , 

4. Description/Purpose: (5. Method; 6. Restriction/ Range; When Applicable) Computes the 
equilibrium conditions for the combustion of any hydrocarbon in air or in any oxygen- 
nitrogen mixture of a given ratio. Uses method devised by Huff, Gordon and Morrell (1) 
involving solution of equations simultaneously by successive iterations. Restricted to 
adiabatic combustion processes in which the total number of elements ' and compounds does 

not exceed 7 and 1% respectively. Thermodynami c functions Program is restricted to 

use with the Table of Coefficients included with the writeup. . 



7. Specifications (Check or fill in appropriate spaces): 

a. Storage used by program: 20K 

b. Equipment required by program: 

Card System x » Magnetic Tape System ; No. of Tapes 



Paper Tape System s Disk File System ; No. of Packs 



TNS, TNF, MF ; Auto divide ; Indirect addressing ;Floating point hardware 



Other (specify) 



Can program be used on lesser Machine ? jj n Specify which requirements can be 

easily removed None 

c. Programming type (Check appropriate spaces): 

Fortran without Format ; Fortran with Format ; 

PDQ FORTRAN with Format X_; 
Fortran II ;/Mainline, Complete X ; Subroutine or function subprogram(S or F) ; 

Is the program a library (ie, SPS) function to the Fortran system checked? He- « 

SPS ; SPS - 1620/1710 ; 



Mainline, Complete ; Macro ; Subroutine 



Other programming language: ; Give details 



d. Language used in the writeup: English 

8. Additional Remarks: The program consists of three parts referred to in the writeup as 
Thermodynamic Functions Program. Combustion Process I and Combustion Process II. 
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Deck Key 



Deck 1 



FORTRAN Source Deck of Thermodynamic 
Functions Program 



Deck 2 



FORTRAN Source Deck of Combustion 
Process I Program 



Deck 3 



FORTRAN Source Deck of Combustion 
Process II Program 



Deck U 



Compressed Object Deck of Thermodynamic 
Functions Program Compiled by PDQ 
FORTRAN Processor C2 with Fixed Format 
Subroutine 



Deck 5 



Compressed Object Deck of Combustion 
Processor I Program Compiled by PDQ 
FORTRAN Processor C2 with Fixed Format 
Subroutine 



Deck 6 



Compressed Object Deck of Combustion 
Processor II Program Compiled by PDQ 
FORTRAN Processor C2 with Fixed Format 
Subroutine* 



Deck 7 
Deck 8 



Sample Problem No. 1 (Input) 



Control Card for Sample Problem Ho. 1 
(One-Card-Input Data for Combustion 
Process II) 



Deck 9 



Thermodynamic Coefficients for Temperature 
Range of l500°-2500° for the Following 
Species: Hg, HgO, OH, C02, CO, H2, HO, 
°2> H » c > H > AR» Used for Sample Problem 
Ho. 1 



Deck 10 Thermodynamic Coefficients for Temperature 

Range of 2$00 o -3500° for the Following 
Species: H 2 , H 2 0, OH, CO2, CO, N 2 , HO, 
O2, H, C, H, AR, Used for Sample Problem 
Ho. 2 

* - First two (2) cards which seroes core hare been removed. 



1.0.0 



e 
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Deck 11 
Deck 12 

Deck 13 



Sample Problem Mo. 2 (Input) 



Control Card for Saaple Problem Ho. 2 
(One -Card-Input Data for Combustion 
Process II) 



Thermodynamic Coefficients for Temperature 
Range of l$00°-2$00° for the Following 
Speciesr H 2 , H 2 0, OH, C0 2 , CO, N 2 , NO, 
O2, HoS, S 2 , SOg, S0 3 , H, C, N, AR, 0, S 
Used for Sample Problem No. 2 
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PROGRAM DESCRIPTIONS 



2.0.0 



COMPUTATION OF EQUILIBRIUM COMPOSITION AND TEMPERATURE 
OF CHEMICAL REACTIONS 



GENERAL DESCRIPTION 



The 1620 Program computes the equilibrium conditions for the 
combustion of any hydrocarbon in air or in air/ oxygen-nitrogen mixture 
of a giTen ratio. The final solution is expressed in terms of the 
partial pressures of each of the reaction products, "P^" , the equi- 
librium temperature "T", and the number of formula weights of the 
reactants, "A" , which are involved in the chemical reaction. 

A rapidly convergent approximation process that simultane- 
ously determines both composition and temperature resulting from a 
chemical reaction is the method used in the program. This method was 
developed by the NACA Lewis Laboratory during 19hB and is applicable to 
a wide variety of problems. 

The program, written for the 1620 computer in IDRTEAN language, 
consists of three passes: (l) Combustion Process I, (2) Combustion 
Process II and (3) Thermodynamic Functions. Input and output data are 
on cards. However, the final solution to a problem will be typed out 
on the 1620 console typewriter. 



© 



2.1.1 



THERMODYNAMIC FUNCTIONS 



PROGRAM DESCRIPTION 



The program calculates thermodynamic data at any temperature 
T (ranging from 300°-3500° K) for use in the Combustion Process I 
Program. 

Thermodynamic data obtained from the program will be the 

following: 

(1) C° - specific heat at constant pressure and standard 

conditions. cal./(raole, °K) 

(2) HfJ - sum of sensible enthalpy and chemical energy at 

temperature T and standard conditions. (Kcal/mole) 

(3) AH - enthalpy change divided by gas constants times 
RT temperature multiplied by -1. 

(U) Log K - logarithm of equilibrium constant. 

(5) S° * - molar entropy at standard conditions. cal./(mole, °K) 



o 

* S T is included in the output but is not used in the Combustion Program. 

Input data to the program are the Coefficients of Thermo- 
dynamic Tables II, III, and IV. Further discussion of the tables will 
be given under the sections entitled "Mathematics" and "Restrictions". 

All input and output data in the program will be on cards. 



2.2.1 
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MATHEMATICS AND/OR TECHNIQUES 



3.0.0 




MATHEMATICS AMD/OR TECHNIQUE 



The method used in the program for the computation of equi- 
librium temperatures and compositions of adiabatic combustion reactions 
was devised by Huff, Gordon and Morrell.(l) A set of equations, repre- 
senting dissociative equilibria, pressure balance, mass balance and 
energy balance is solved simultaneously by successive iterations of a 
series of correction equation. Each of these correction equations con- 
tains a correction variable which takes the form: 

&x « Ax log £o (1) 
X 

Due to the limited storage capacity of the 1620 computer and 
the necessity to calculate thermodynamic data for each iteration, the 
program is written in three parts. The function of each part is as 
follows: 

(1) Using the estimates of the partial pressures of 
the reactants involved in the combustion process, 
thermodynamic properties determined by the Thermo- 
dynamic Functions Program, an estimated temperature, 
the general formula and total enthalpy as basic 
data, Combustion Process I (Deck 1) constructs 

two matrices in the manner described by Huff, 
Gordon and Morrell,(l) and tests for convergence 
using .0015 as the tolerance level* 

(2) Combustion Process II (Deck 2) multiplies the two 
matrices that were constructed in Combustion 
Process I and computes new estimates from the 
matrix obtained as a result of the matrix multi- 
plication by using the Croat* ^) method for 
solving a matrix , as suggested by Huff, Gordon 
and Morrell.(l) 

(3) The Thermodynamic Functions Program (Deck 3) 
computes the thermodynamic properties given in 
equations (2) - (6) for use in the Combustion 
Process I Program (Seek 1). The coefficients 
and constants in Tables II, III and IV used in 
computing the thernodyn&mic properties necessary - 
for Combustion Process I were computed specifi- 
cally for the purpose of being used with the 
Thermodynamic Functions Program. The following 
equations representing the molar specific heat, 
the molar enthalpies, the molar entropies, the 
enthalpy change and the equilibrium constant 

are used in the program: 



3.1.1 



e 



C° '<=<* flT +Tt 2 + «T 3 (2) 

p 

H o » (e< T + i/2gT 2 + 1/3 "if T 3 + 1A6T 1 *) 10-3 + C (3) 

T ^ 

S° ■= °< lrff + PT + 1/2TT 2 + 1/3ST 3 + K (U) 

AH - a + bT + cT 2 + dT 3 + m (5) 

HF T X 10-3 

Log K « A log T ♦ BT +T T 2 + &T 3 + L + M (6) 

T X 10-3 



Data contained in the National Aeronautics and Space Adminis- 
tration Report (3) were used to determine the coefficients in equations 
(2) - (6). However, applying equation (l) to the conservation of energy: 

\ - Ah log ho (7) 
h 



where "h " is the enthalpy of the reactants and "h"- is the sum of the 
enthalpies of the product weighted by their respective mole fraction of 
the combustion products. If for any reason, the values of "h " and "h" 
would be opposite in sign, the expression containing the log (h /h) would 
create an error response in the computer. Therefore, to evaluate "C" of 
equation (3), assigned enthalpies, as shown in Table I, were used. 



3.1.2 
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TABLE I 

ENTHALPY ASSIGNED TO SEVERAL SUBSTANCES 



Element 


Phase 


H298 


Ar 


Gas 


1.U812 


Br 2 


Liquid 


2.6600 


C 


Graphite 


92.1790 


H 2 


Gas 


69.UU07 


N 2 


Gas 


3.7715 


°2 


Gas 


U.1109 



S Rhombic 100.0000 



3.1.3 
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INPUT AND OUTPUT FORMATS 
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INPUT /OUTPUT FORMAT 
THERMODYNAMIC FUNCTIONS 



a 



INPUT/OUTPUT FORMATS 



THERMODYNAMIC FUNCTIONS 



Input 

A set of input data for a compound or element being processed 
through the Thermodynamic Functions Program consists of three cards 
(one from each of the Coefficients of Thermodynamic Functions Tables). 
A 2, 3 or h will be punched in Column 80 for Tables II, III or IV, 
respectively. A set of data being processed wast be in numeric order 
by table code. 

All data being processed are preceded by a temperature card. 
The temperature card will be the first card read into the program and 
will have the following format: 



Description Columns FORTRAN Symbol 

T (°K) 1-8 T 

-503059 069916 -006240 000319 385222 4775775 
503059 -160988 028738 -002206 10996629 
495233 638682 -172814 018124 1230309 1580923 
-207869 016067 -000377 -000032 090789 2246215 
207869 -036996 001736 000218 5172101 
580531 147043 -010349 -001733 6768494 -180405 
2200.00 TEMPERATURE CARD 



Format 
F8.2 

H20 
H20 
H20 
H2 

H2 
H2 



j j , , a , o c o o o o o a o o o c o , o o o o o o , o o o o o o o o o o o o o o o o o o o o n o o o o o o o o o o o b o o o o o o o o o o o o o o o a a 

1 j d 7 S 3 ;0 -I II :'l M : 5 IS C :h 'S ?C I" li::2* T-'~ I". 7B 29 3C y. T ;:6 ?7 38 39 40 4142 4344 454S47 W 49 50 51 52 53 54 55 56 57 H 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 7* 75 76 77 78 79 30 
11 1 1 1 1 11 il 11 11 MMM 1 i 1 1 i , : , • \ \ lo f*WWkl 111111 1111111111111111111111111111 





2 , i 2 2 2 2 2 2 2 2 2 2 2 !»Wf\ 2 2 fife 2 2 
33333^33333333 *****T3 i /3 3 3 l3 A »3 3 3 

44444444444444444444444k f**<4 & 4 A4 4' 

5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 E 5 L 5 5| 5j^\s 5 5 I 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 

6 6 6 6 6 6 6 6 ( 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 Mfi/i 6 6 16 6 1 6 6«f 6 6 6»6 S 6 6 6 6 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 



2|2 2 2 2 2 |M2 22222222222222222222222222222 
3 3 3 3 3 MM**3 ftftM*3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 

4 A* 4444444444444444444444444444444444444 



77777777777 777777 7777777,7777777777777777777777777777777777777777777777777 777777 
8 8 8 8 8 „ 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 » 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 6 8 8 I 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 
99999999999999999999999999_999_9999_99999999999999999999999999999999999999999999 

1 2 3 4 5 6 7 t 9 10 I I 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 21 » 30 il 32 33 34 35 36 37 31 39 40 41 42 43 44 45 46 47 48 49 SO 51 52 53 54 55 56 57 58 59 fit) 61 12 63 S4 65 66 67 68 69 70 71 72 73 74 75 76 77 71 79 80 
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Input data from all three tables are on cards and will be in 
the following formats: 

Table II 



Description 
I x 103 



106 
10? 



Identification 
Table Code 



Columns 

1-8 
9-16 
17-21 
25-32 
33-4;l 
U2-J*9 
55-57 
80 



FORTRAN Symbol 

A 

B 

CC 

D 

E 

F 

ID 



Fontat 

F8.5 
F8.5 
F8.5 
F8.5 
F9.5 
F8.5 
12 



Table III 



Description 



b x 



10 i 

c x 10° 
d x 10? 
m 

Table Code 



Columns 

1-8 

9-16 
17-2U 
25-32 
33-U1 
80 



FORTRAN Symbol 



AA 

BB 

CCC 

DD 



Pd mat 

F8.5 
F8.5 

F8.5 
F8.5 

F9.5 



Table IV 



Description 

B x^O* 
C x 1 6 
ft x 10? 

L 

M 

Table Code 



Co Irani 

1-8 

9-16 
17-21; 
25-32 
334a 
U2-l*9 
80 



FOBTBAH Symbol 



CK 
DDD 



Format 

F8.5 
F8.5 
F8.5 
F8.5 
F9.5 
F8.5 



U.1.2 
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Output 



All output data are on cards. There will be one card punched 
for each set of data read. The order of the elements and compounds in 
the output is the same as their original order in the input. 



The output format will be as follows: 
Description Columns FORTRAN Symbol 



cal./(raole, °K) 
(leal/mole) 



1-9 



10-18 



19-27 



HH 



Format 



F9.U 



F9.U 



F9.lt 



log K 

cal./(raole, °K) 
Identification 



28-36 

37-U5 
76-80 



CBCF 

S 
ID 



F9.U 

F9Jt 
15 
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148721 


-061149 


017015 


-001596 


5108413 




S2 


3 










SAMPLE INPUT DATA 
FOR 








-148721 
111 0546 


026557 
292319 


-003695 
-103367 


000230 
013607 


-156036 


2218556 

-D133 U^: 


S2 
S02 


4 

c. 




- 




THERMODYNAMIC FUNCTIONS PROGRAM 








33444 


- 098 02 


0238 69 


-0022 09 


1 *5 O O "2 O O /"\ 




S02 


3 


2200.00 








TEMP. = 2200 K 








-233444 


042570 


-005183 


000320 


-642 058 


5595112 


S02 


4 


58 0531 


147043 


-01 0349 


-001733 


6768494 -180405 




H2 


2 




A A 1*1? 
H-Dut DO 


— 1 71 2 84 


02 2 5 66 


58 9454- 


^) UO ZJ £. O 


S03 


2 


207869 


-036996 


001 736 


000218 


5172101 




H2 


3 


2 779 04 


— 1 4 05 7 


035 01 2 


-UUJt D X 






• S03 


3 


-207869 


016067 


-000377 


-000032 


090789 2246215 




H2 


4 


-277904 


061 049 


-007603 


000472 


-1343033 


7417350 


S03 


4 


4952 33 


£ o o c o o 


— 1 728 14 


018124 


1230309 1580923 




H20 


2 


aq f^p i n 


noonno 

U \J \J \J \J \J 


uuuuuu 


nnnnnn 

U \J\J \J \J \J 


O _J 3 3 D ^2 U 


- 091 377 


(_l 


2 


5 03 05 9 


-160988 


0c. o I io 


0022 06 


10996629 




H20 


3 
















3 


-503059 


069916 


-006240 


000319 


OODd^C. 4 / I J 1 ID 




H20 


4 














H 


4 




o y, a c o q 


— 061 041 


005829 


4467687 1312034 




OH 


2 


515 5 97 


— 02 8 602 


01229 






QCIOOO A 
O 30 O C H- 


Q 


2 


2 31 047 


-0615 66 


009990 


-000659 


5073363 




OH 


3 














Q 


3 


-231047 


026738 


-002169 


000095 


188357 2203334 




OH 


4 














C 


4 




5 02 247 




022673 


141059 -643266 




C02 


2 


489391 


01322 6 


-0078 69 


0015 68 


1 1346040 


872343 


N 


2 


279290 


-1 34150 


0312 09 


-002810 


19305350 






3 














N 


3 


-279290 


058261 


-006777 


000407 


-639973 8384207 




C02 


4 














N 


4 


620126 


255577 


-088801 


011249 


6578319 1103990 




CO 


2 


49681 


uuuuuu 


UUUUUU 


UUUUUU 


/~i i-i n n An 
UUUUUU 




A D 

AK 


2 


199668 


-071794 


016708 


-001447 


12896790 




CO 


3 












/. .- 


AR 


3 


-199668 


031180 


-003628 


000210 


-074022 5601005 




CO 


4 














AR 


4 


595705 


271090 


-092791 


011599. 


186228 1106826 




N2 


c. 


5 01322 


-001 167 


— 001484 


000591 


601608 9 


1 000671 





2 


192766 


-061552 


012924 


-001065 


11325067 




Vic. 


3 

















3 


-192766 


026732 


-002806 


000154 


-055680 4918414 




N2 


4 

















4 


664779 


221014 


— 077591 


009955 


2332637 1159108 




NO 


2 


5 71818 


— U^*4- 3 u 




nnc i ^ n 


1 OH-O J. / DU 


787223 


c 


2 


1 6401 2 


— 052 574 


011446 


-000981 


75 62922 




NO 






- 










s 


3 


-164012 


022832 


-002485 


000142 


-084787 3284535 




NO 


4 














s 


4 


7612 06 


no/, a n7 


— 01 6012 


002108 


148991 472934 




02 


2 


















121493 


-024340 


002188 


-000117 


5979684 




02 


3 


















- 12 1493 


010571 


— UUU4 ( i> 


000017 


-281254 2596944 
























o06270 


702346 


-230861 


028578 


16243440 1252268 




H2S 


2 


















482659 


-200597 


045752 


-004241 


8708322 




H25 


3 


















-482659 

848 089 


087118 
053176 


-009935 
-017635 


000614 
002345 


351185 3781976 
22811776 564902 




H2S 

S2 l?S 


4 
2 
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INPUT/OUTPUT FORMAT 
GOMBDSTIOH PROCESS I 



INPUT /OUTPUT FORMATS 



COMBUSTION PROCESS I 



Input 

A set of Input data into Combustion Process I will be in the 
following format and order: 



1. Control Data Card (One Card) 
Card Column 



1-3 Number of Compounds 

U-5 Number of Elements 

78-60 101 (Identification No.) 



Format 

13 
12 



FORTRAN Symbol 



2. A, P, T Data Card (One Card) 
1-11» 
15-28 



29-U2 
76-80 



Number of Formula Weights Eli*. 7 

of Reactamts 
Initial Estimate of Total ElU.7 

Pressure 
Temperature ELU.7 
201 (Identification No.) 



3. Log ^(i) Data Cards (No. of Cards 
by No. of Reactants) 



rained 



l-lli Log P(i) 

78-80 301, 302, Etc. 

(Identification No.) 



ELU.6 



B(l) 
P(D 
T(l) 



K8,i) 



U. Thermodynamic Data Cards (No. of Cards Determined 
by So. of Reactants J 



1-9 


o 

C 

P 


F9Ji 




10-18 


H° 
T 


F9.U 


1(9,1) 


19-27 


iH/RT 


F9.lt 


K9,D 


28-36 


Log K ± 


f9.h 


1(10,1) 


37-U5 


S° 
T 


F9.U 




78-80 


U01, U02, Etc. 







U.2.1 



o 



5. a , b c , , hj, Data Card (One Card) 



Card Column 
1-9 

10-18 

19-27 

28-36 

37-U5 

U6-5U 

55-63 

6U-72 

78-80 
6* a^, b^, • • 



a Q - Mo. of Atoms of 
Respective 
Element in 
Formula 

b c - No. of Atoms of 
Respective 
Element in 
Formula 

c c - Ho. of Atoms of 
Respective 
Element in 
Formula 

dj, - No. of Atoms of 
Respective 
Element in 
Formula 

e Q - No. of Atoms of 
Respective 
Element in 
Formula 



f„ - 



8o " 



No. of Atoms of 
Respective 
Element in 
Formula 

No. of Atoms of 
Respective 
Element in 
Formula 

Initial Enthalpy 
per Mole of 
Reactants 



Format 
F9.lt 



F9.lt 



F9.lt 



F9.lt 



F9.U 



F9.U 



F9.lt 



F9.lt 



FORTRAN Symbol 
0(1) 

0(2) 

0(3) 

0(10 



^01 (Identification No.) 
,, Data Card (No. of Cards Determined 



by No. of Reactants) 

1-9 a^ - No. of Atoms per 

Element in 
React ant Product 

10-18 bi - No. of Atoms per 

Element in 
Reactant Product 



F9.U 
F9.lt 



0(5) 



0(6) 



0(7) 



0(6) 



T(l,i) 
T(2,i) 



U.2.2 



o 



6. b ± , . . ., gi Data Card (No. of Cards Determined 

by No. of Reactants) (Cont.) 

Card Column Format FORTRAN Symbol 

19-27 c^ - No. of Atoms per F9.U T(3»i) 

Element in 
Reactant Product 

28-36 d^ - No. of Atoms per F9.lt T(U»i) 

Element in 
Reactant Product 

37-U5 e ± - No. of Atoms per F9.U *(5»i) 

Element in 
Reactant Product 

U6-5U f • - No. of Atoms per F9.lt *(6,i) 

1 Element in 

Reactant Product 

55-63 g± - No. of Atoms per F9.it T(7,i) 

Element in 
Reactant Product 

78-80 601, 602, Etc. (Identi- Q(i) 

fication No.) 



U.2.3 



Output 

Output data for Combustion Process I will be from the 
typewriter if the error parameter "E" is less than or equal to .0015* 
If E is greater than .0015* the data in core store are left in the 
machine and are used as data for Combustion Process II. The typewriter 
will type out the following message: 

"Push reset, Load Combust* Prog* 2." 

A sample output listing of Combustion Process I is shown on 
the next page. 



© 



U.2.U 



SAMPLE OUTPUT 



COMBUSTION PRO CESS I 





H2 


. 11437648E-02 




H20 


.73096864E-01 





OH 


.23769886E-02 





C02 


.15 166840 E 00 





CO 


.12786094E-01 


" - — 


N2 


.73562902E 00 




NO 


.29623556E-02 


. 


02 


.9211 4445 E-0 2 




H2S 


. 17641 303 E-#8 




S2 


. 20 1 87753 E- 10 




S02 


.17414738E-02 




S03 


.50231170E-06 . j 




H 


.24336940E-03 / 




C 


.32717617E-16 




N 


.17737053E-07 




AR 


.8801 926 1E-02 







.3406953 1E-03 


20 
19 


S 


* 58250858 E-07 


18 

17 


16 
15 


A= 


.30561900E 00 


14 
13 


T= 


. 22448700 E 04 


12 
11 


STOP 0000 



10 
9 
8 
7 
6 
B 

a 

3 
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INPUT /OUTPUT FORMAT 
COMBUSTION PROCESS II 







INPUT /OUTPUT FORMAT 
COMBUSTION PROCESS II 



Combustion Process II requires only one data card. The data 
card is a control card and has the following formats 

Card Colmm Format 

1-3 No. of Compounds 13 

No. of Elements 12 

Output 

Output data will be on cards and will hare the following 

formats 

Card Column Format FORTRAN Symbol 

1. Temperature Card 

1-8 F8.2 T ( 1 ) 

2. Control Data Card 

1-3 No. of Compounds 13 * 

U-5 No. of Elements , 12 M 

78-80 101 (Identification/ 

No.) / 

3. A. P, T Data Card (One Card) 

1-llt No. of Formula Weights EUt.7 B(l) 

of Reactants 

15-28 Initial Estimate of ElU.7 P(l) 

Total Pressure 

29-U2 Temperature Elli.7 T(l) 

78-80 201 (Identification 

No.) 

U. Log P(i) Data Cards (No. of Cards Determined 
by No. of Reactants) 

1-llt Log P(i) K8,i) 

78-80 301, 302, Etc. (Identi- 

fication No.) 



U.3.1 



RESTRICTIONS 



5.0.0 



RESTRICTIONS 



The program is restricted to computation for adiabatic chemical 
reactions in which the total number of compounds and elements does not 
exceed 13 and 7, respectively. A code card will allow any variable 
number of compounds and elements to be processed. (For example, 5 
compounds and 2 elements, 11 compounds and 5 elements, 12 compounds and 
7 elements, etc.) The code card is also used to avoid division by zero 
and to avoid the possibility of taking the log of a zero quantity. 

In the test for convergence, the program requires that e * 
.0015. However, the smaller the value of e, the more accurate the final 
results will be. Thus, for a greater degree of accuracy the user of the 
program may want to require e to be smaller than .0015. 

The coefficients and constants in Tables II, III and IV used 
to determine thermodynamic properties are limited to a total number of 
seven elements and seventeen molecular or radical species. If thermo- 
dynamic data for other species are required, the coefficients and 
constants in equations (1) - (5) may be determined by referring to 
reference (3). However, to determine "C in equation (2), reference 
should always be made to the assigned enthalpies in Table I. If an 
element not listed in Table I should be involved, a value may be as- 
signed to represent its molar enthalpy in its standard state. The 
assigned value should be selected so as to result only in positive 
values of enthalpy for all other species containing this particular 
element* 

The constants and coefficients in Tables II, III and IV are 
for use only with the Thermodynamic functions Program submitted with 
this writeup. The thermodynamic values obtained from the program by 
the use of these equations are not internally consistent with values 
obtained from other sources. 

If the original estimates are too unreasonable, the results 
obtained in any given iteration may exceed, in temperature, the thermo- 
dynamic equations furnished with this program. In such cases, it is 
suggested the program be reconsidered with more reasonable estimates. 



* - The FORTRAN symbol for e is E. 



5.1.1 
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OPERATING INSTRUCTIONS 



OPERATING INSTRUCTIONS 



THERMODYNAMIC FONCTIOtE 



Sorting Input 

Input data are in the order as indicated in the Input/Output 
Formats for Thermodynamic Functions Program. No further sorting is 
required. 

Clearing Core 

The Object Deck, which was compiled by FDQ FORTRAN, auto- 
matically clears core when it is loaded in the machine. 

Switch Settings 

All program switches are set to "off. 

Loading Object Deck and Data 

Place Object Deck in 1622 read hopper, followed by tempera- 
ture card and data punched from Thermodynamics Coefficients Table.* 
Push "load" button. 

Program Stops 

There are no "halts" in the program. 

Output Data ^ 

All output data are on cards and require no sorting. After 
all input data have been processed, the machine's "reader no feed" 
light will be on. Retain output data for input into Combustion Process 
I Program. 



*Note: The set of data used depends on the temperature. Each time a 
new temperature is used, the program must be reloaded, followed 
by the temperature card and data from the tables for that 
particular temperature range. 



6.1.1 



OPERATING INSTRUCTIONS 



COMBUSTION PROCESS I 



Sorting Input 

Input data are in the order as indicated in the Input/Output 
Formats for Combustion Process Program I. No further sorting is re- 
quired. 



Clearing Core 

The Object Deck, which was compiled by PDQ FORTRAN, auto- 
matically clears core when it is loaded in the machine* 



Switch Settings 

All program switches are set to "off". 



Operation (1st Iteration) 



1. All data for the first run should be carefully 
edited. 



2. Load Object Deck. 



3. Load data* 



It. If the following message is typed out on the 
typewriter: 

"Push reset, Load Combust. Prog. 2." 

Do not clear the machine. Push reset and load 
Combustion Process II. If the desired answers 
to the problem have been obtained, the final 
answers will be typed out on the typewriter in 
the format indicated in the writeup under Input/ 
Output Formats for Combustion Process I. 
Processing will then be complete. 

Operation (2nd Iteration) 

1. Remove first four sets of original input data. 



2. Replace original data with new data obtained 
from Combustion Process II, followed by new 
data obtained from Thermodynamic Functions 
Program. 

3. Repeat Steps 2 - h of first iteration. 
Operation (3rd, Uth, $th, Etc. Iterations) 

1. "Same as second iteration". 

6.2.1 




Program Stops 

The program will come to a halt only when processing is 

complete* 



6.2.2 

o 




OPERATING INSTRUCTIONS 
COMBUSTION PROCESS II 

Input Data 

Input data are left in machine from Combustion Process I. 
One data card ( control data card) is also read in as input to the 
program* 

" Do not clear machine before loading Program" . 

Loading Program 

When message is typed out from Combustion Process I, the 
program will come to a halt. Push reset and load program with data 
control input card on back. 

Output 

All output data are on cards. Retain temperature card (1st 
output data card) for processing through Thermodynamic Functions Program. 
The remaining data will be processed through Combustion Process I. 

Program Stops 

The program comes to a halt only when processing is complete. 



6.3.1 



a 



a 



SUPPORTING DOCUMENTATION AND MATERIALS 



7.0.0 



SAMPLE PROBLEMS 



SAMPLE PROBLEM NO. 3 
INPUT DATA 



8 5 

0.0000001E+07 0.00O0QO1E+Q7 . 00022 00E+07 



-3 . OOOOOOE+OO K2 
-2 .732828E+0O H20 
-2 .698970E-! 00 OH 
-2 .031517E+00 C02 
-2.698970E+0Q CO 
-2.148742E+00 N2 
-3. OOOOOOE+OO NO 
-3 . OOOOOOE+OO 02 
-4. OOOOOOE+OO H 
-15 .OOOOOOE+OO C 
-6, OOOOOOE+OO N 
-2 . 09691 OE+00 AR 
-4. OOOOOOE+OO o 

8.3548 83.5462 -24.8815 4.5019 45.7977 
12.5689 33.5819 -53.0993 10.0160 64.4355 

8.4137 60.6026 -24.4301 4.6695 58.6888 
14,5416 29.8104 -88.8044 23.3723 75.2756 
8.7237 83.1178 -59.6935 18.5778 62.6390 

8.6650 18.9139 -52.5632 15.8252 61.0450 

8.8147 41.1290 -35.3099 8.9970 66.0917 

9.1384 20.0762 -27.9532 5.1^3 65.0775 

4.9681 96.2667 -.0000 .0000 37.3217 

5. C315 272.9527 -.0000 .0000 47.7078 

4.97 09 12 4.3596 -.0000 .0000 46.5441 
±^966 1_ 10_.9298 -.0000 .0000 46.9124 

4.9786 71.1436 -.0000 .0000 48.5488 7.2.3 



o a 

4. 0000 1 .00? Q 14.9090 0.0890 4^0057 — - 24 9. 6473 

-2U.J_000 ... _.: _ _ H2 

2- 000 1.0000 . . ... H?n 

1-0000 .... 1.0000 . OH 

1. .00.00 : 2.0000 C02 

L_ QQ0O ; ....... ■ _ . 1.0000 ... ... CO ... 

: : 2^0000 • H2 

1___00 1.0000 NO 

- . .... 2-.0QOO . .. ._£_i .... 

1.0000 jd 

l . oooo c 

1.QQQ0 __ 

. 1.0000 . Ag 

1 .0000 . O 



17 



7.1.2 



a 



SAMPLE PROBLEM NO. 1 
ANSWERS 



H2 1 .3375881 8E-02 

H20 .183^333 IE 00 

OH . 2632.5090 E-02 

CO 2 . ._J60 8526^0 1__ 

CO .8351 1769E-02 

N2 .701£i6i 10 E 00 

NO . 18175524E-02 

02 .42697 982 E-02 

H .34202367E-03. 

C . 1 50 100 1 5E-1 6 

N . 1 132326 IE-07 

AR .83 80 6500 E-02 

.1.85J0828E-03 

A- .941 63000 E-01 

T= .22085000E 0k 

STOP 0000 



./ 
7.1.2 



SAMPLE PROBLEM NO. 2 
INPUT DATA 



12 6 

. 00OO031 E+05 0.QO000O1E+07 . 00022 OOE+07 

-2.52288E+00 

-1 .22185E+O0 

-2.69897E+00 

-0.8 2391 E+00 

-1.69897E+00 __ 

-0.I5490E+00 _ _ 

-2.69897E+00 

-2.301 03E-fOQ 



■6.00000E+00 



-6.00000E+QQ 



-2.69S9I7E+00 



■6.0000QE+00 



-3.39794E+Q0 



-17.0000E+00 



-5.30103E+00 
-2.05061E+00 
-3.69897E+00 



-8.0Q0OQE+00 _ 
8.3548 83.5462 -24.8815 4.5019 45„7977 

12.5689 33.5819 -53.0993 10. 0160 _64.4355 

8.4137 60.60 26 -24.4301 4.6695 58.6888 

14.5416 29.8104 -88.8044 23.3723 75.2756 



.7237 83.1178 -59.6935 18.5778 62.6390 



8.6650 18.9139-52.5632 15.8252 61.0450 



!.8147 41.1290 -35.3099 8.9970 66.0917 7.2.1 



_9^13B.4 2Q.0Z£2^ 



13.3836 
9,0469 
_13.9823_ 

19^4736_ 
4 , 968J. 
5. 0315 
4.9709 
4.968 1_ 
4.9786 
5.1115 
0.4944 



186.2487 
247,5739 
57.7504 
45.7924 
_96.2 667 
272.9527 
124.3596 
10.9298 
71.1436 

0.5381 



2.0000 



2.0000 



-27^9532 


.. 5 


.1423 


.65^0775 


-41.7595 


6 


.0712 


70.0615 


-24,0155 


3 


.9830 


71.7460 


-59,6581 


11.9288 


83.7493 


-78 ,6660 


12 


.0215 


94,9493 


-.0000 




.0000 


37.3217 


OQOO 




. OOQO 


47,7078 


- . 0000 




.0000 


46,5441 


-.0000 




.0000 


46.9124 


-,Q000 




. 0000 _ 


48.5488 


-.0000 




.0000 


50.6240 


4,8237 





,0288 


1,.?638 








1^0000 



1.0000 



I. 0000 



1.0000 



2..0000 
1. 0000 



1.0000 

__2.0O00. 
1.0000 

1 . 0000 
2.0000 



2 .0000 



0,0057 



1 . 0000 



2.0000 
2.0000 1.0000 
3. OOQO 1.0000 



1 . 0000 



1 . 0000 



1 . 0000 



1.0000 



1 . 0000 



1.0000 



75^8600 

H? 

H20 

OH 

C02 

CO 

H2 

02 

_.H2S 

S2 

S02 

S03 

H 

C 

N 

AR 



S 



7.2.2 



© 



a 



SAMPLE PROBLEM NO. 2 
ANSWERS 



H2 . 1 1437JS48E-02 

H20 .73096364E-01 

OH .2376^386E-02 

CO 2 .15 166340 E 00 

CO . 12/36034E-01 

N2 .73562902E 00 

NO . 29623 556E-02 

02 .9211 4445 E-02 

H2S . 17641 303 E-08 

S2 . 20 1 8775 3 E- US 

S02 .174 147 38 E-02 

SQ3 .502311 70 E-06 

H .2433694 0E-03 

C .32 717617E- 16 

N . 1773 7053 E-07 

AR _ .880 19261 E-02 

'l k ®^ 9 l 3 J3^ 

S • 58z i ( ^58E"?_Z 

A= .30561 30$ E 00 
T= . 22448700 E 04 
STOP 0000 



7.2.3 



SOURCE PROGRAM LISTINGS 



7.3.0 



PROGRAM LISTING 

THERMO DYNAMIC FUNCTIONS 



2600008000(995 
PDQ FORTRAN C2 
START 



-6600 READ 30, T _ 

-6624 5 READ 10, A,B,CC,D,E,F, ID _ 

-6720 B=B/10.**3'1 _ 

-6768 CC=CC/10.**6. _ 

-6816 D=D/t0'.**9. . _ 

-6864 X=0._43_4294*LOG(T) 

-6900 C=A+B*T+CC*T**2'.+D*T**3?. _ 

-7068 H= ( ( A*T+. 5*B*T**2 .+CC*T**3 ./ 3 .+. 25*D*T**4 J / 1 t^**3_ . )+E 

-7428 S=2 . 302585*A* X+B* T+. 5*C C*T**2 »+6*T**3 . / 3 .+F * 

-7704 * READ 10, AA,BB,CCC,OD,EE " " " 

-7776 BB=BB/10.**3. 

-7824 CCC=CCC/10.**6. 

-7872 DD=DD/10.**9. " _ 

-7920 HH=AA+BB*T+CCC*T**2.+DD*T**3 .+ (EE/ (T* 10.^-3. )) 

-8160 HH= (^-1*HH __~ " 

-8208 READ 10, AAA,BBB,CK ? DDD,EEE,FF _ 

-8292 BBB=BBB/10.**3. 

-8340 CK=CK/t0-.**6. " 

-8388 DDD=DDD/10v**9. 

-8436 CECF=AAA*X+BBB*T+CK*T**2.+DDD*T**3.+EEE+(FF/(T*10.**-3.)) 

-8724 ID=ID+400 " " '_ 

-8760 PUNCH 20, C,H,HH,CECF,S.,ID 

-8844 GO TO 5 

-8852 10 FORMAT (4F8.5,F9.5 s F8.5,5X,t2) 7.3,1 



-8912 20 FORMAT (5F9.4,30X, 15) 
-8966 30 FORMAT (F8.2} 
-8988 END 
T9999 SIN 
T9989 SINF 
T9979 COS 
T99&9 COSF 
79959 EXP 
T9949 EXPF 
T9939 LOG 
T9929 LOGF 
79919 SQRT 
T9909 SQRTF 
T9899 ABS 
T9889 ABSF 
T9879 DRH 
T9869 DRHF 
T9859 ATAN 
79849 ATANF 
T9839 0030 
79829 T 
79819 0005 
79809 0010 
79799 A 
79789 B 
79779 CC 
79769 D 
79759 E 

79749 F 

79739 ID 
79729 521 

79719 51 

79709 000 7.3.2 



19699 51 < 
T9689 511 

T9679 X 

T9669 5043429400 
T9659 C 

79649 5120000000 
79639 301 
79629 002 
79619 H 

79609 5050000000 
79599 003 
79589 004 
79579 005 
79569 5025000000 
79559 006 



79549 51400 00000 
79539 007 

79529 008 

79519 S 

79509 5123025850 

T9!*99 AA 

79489 BB _ 
79479 CCC 
79469 DD 



79459 EE 
79449 HH 

7943? 51; 

79429 AAA 

79419 BBB 

79409 CK 

79399 DDD 

79389 EEE _ 

79379 FF 7.3.3 



a 



79369 CLCF 
79359 0400 
79349 0020 



LOAD SUBROUTINES 
260000800009RS 



PROGRAM LISTING 
COMBUSTION PROCESS I 



PDQ FORTRAN C2 

START 

-6600 C COMBUSTION PROCESS 

-6600 C PROGRAMMED FOR RESEARCH DEPARTMENT 

-6600 C OCTOBER, 1964 _ 

-6600 D I MENS ION X(9,23) , Y( 10 , 23 ) , B( 1 >, P( 1 ) ,T( 1 ) ,0(8) ,Z(9) ,Q (20) 

READ 60,N,M 



-6636 _ M=M+13 

-6672 DO 99 1=1,9 

-6684 Z(t)=0. 

-6720 J_° JL 9 _ _ll* 2 __ 
-6732 99 X(l,J)=0. 

-6876 DO 30 #=1 ,-t0_ 

-6388 DO 30 J=t,23 



-6900 30 Y(f , J)=0 . _ 

-7044 READ 402,B(O,P(lf,TfO 



-7092 DO 400 1=1, N 



-7104 400 READ 402, X (8, I ) 
DO 403 1=14, M 



-7224 403 READ 402 , X (8 , t ) 



-7332 K=1 

-7344 L=N_ _ 

-7356 8 DO 9 hKJ. 

-7368 9 READ 100,X(1 , I ) ,X (9, I ) , Y(9, I ) , Y( 10, I ) 

-7656 IF (L-N) 10,11,10 

-7724 11 K=14 _ 

-7736 L=M 

-7748 GO TO 8 



-7756 10 READ 100,0(1), 0(2>, 0(3), 0(4), 0(5) ,0(6) ,0 (7) , 0(8) 7.3.5 



-7864 DO 404 l=1,N 

-7876 404 X(8,I)=EXP((X(8, I)/. 434294)) 

-8056 DO 405 I=14,M 



-8068 405 X(8,l)=EXP((X(8,|)/.434294)) 

-8248 DO 13 1 = 1,20 ___ 

-8260 X(9,22}=X(9,22)+X(8,t)*X(1,l) 

-8404 X(9,l)=X(8,l)*X(9,l) _ 

-8584 X(9,21)=X(9,21)+X(9,I) 

-8668 13 CONTINUE _ 

-8704 X(9,22)=T(1)*X(9,22)/1000. 

-8752 X(9i2T)=(-1.)*X(9,2T) 

-8800 K=1 . _ 

-8812 L=N _ 

-8824 20 DO 14 l=K,L _ _ __ 

-883 6 1 4 READ 600 , Y ( 1 , 1 > , Y(2,Q,Y(3,0 , Y (4 ,1 )^(5,0^(6,0^ (7,0 ,0.(0 

-9340 IF (L-N) *5 f t7,i5___ 

-9408 17 K=14 ' 

-9420 L=M 



-9432 GO TO 20 

-9440 t5 DO 16 1=1,7 

-9452 DO 16 J=1,20 

-9464 X(r,J)=X(8,J)*Y(l,J) 



-9668 16 X(l,2l)=X(l,21)+X(l,J)_ 
-9884 DO 18 r=l,7 



-9896 18 X(I,21)=X(I,21)*(-1.) > 

J0028 D0_19 1=1,10 

J0040 19 Y(l,l-H3)=1. 

J0148 M=M-13 _ 

J0184 DO 41 J=1,N 

J0196 Y Q J Y (J . j ) + ( LOG ( X ( 8 , J » * . 43 4294 ) 

J0388 DO 41 1 = 1, M __ 

J0400 41 Y(10,J)=Y(10,J)- Y(I,J)*(LOG(X(8, 1+13))*. 434294) 

J0748 C SUM P(l) 7.3.6 



© 



9o 



J0748 
J0760 
J0772 
J0784 
J0868 
J0988 
J 1000 
J 1228 
J 1348 
J 1408 
J 1420 
J 1588 
J 1624 
J 1 684 
J 1696 
J 1732 
J 1744 
J 1756 
J 1852 
J 1920 
J 1932 
J2076 
J2088 
J2232 



J 2 268 
J2280 
J2292 
J2390 
J 240 2 
J2430 
J2462 
J2526 
J2584 



SUMP=0. 
V=0. 

DO 51 1=1,20 



51 



74 



50 



V=V+ABS(Y(10,I)) 
SUMP=SUMP+X(8,I) 

DO 74 l=1_,M 

Z(l)=LOG(B(1)*0(|)/((-1.)*X(1,21)))*.434294 
Z(8)=LOG(B( 1 )*0(8)/_( (-1 . )*X(9^21 )) )*.434294 

Z(9)=L0G(P(J)/SUMP)*. 434294 _ 

DO 50 1 = 1, M 

XO,23)=(-1.)*X(l,20*Z(0^ 

X_(8,23)=SUMP*Z(9) 
X(9,23)=(-1.)*X(9 i 21)*Z{8) 

L=M 

M=M+13_ 

£=V 



83 



DO 83 1=1,9 
E=E+ABS(Z(0) 



200 
210 

310 



212 

60 
402 
600 
100_ 
202 



I F (E-.0015) 200,200,20 1 

D0_ 210 1 = 1 ,N_ 

TYPE 202,Q(I),X(8,|) 



DO 310 l = 1_4,_M 

TYPE 202, Q(l),X(8,i) 

TYPE 205, B(1),T (1) 

STOP_ 

TYPE 212 



FORMAT (/34HPUSH RESET. LOAD COMBUST. PROG . 2. ) 
STOP 

FORMAT (13,12) 

FORMAT C3E14.0) 

FORMAT (7F9.4,9X,A4) 
FORMAT (8F9.4) 



FORMAT (/A4,3H ,E14.8) 



7.3.7 



J2630 205 FORMAT (//2HA=,E14.8/,2HT=,E14.8) 
J2694 END 

T9999 SIN 

79989 SINF 
79979 COS 

79969 COSF 

79959 EXP _ 
79949 EXPF 

79939 LOG _ 
79929 LOGF 
79919 SQRT 

79909 SQRTF _ 

79899 ABS 

79889 ABSF _._ 

79879 DRH 

79869 DRHF 

79859 ATAtJ 

79849 ATANF 

79839 X 77779 

77769 Y 75479 _ 

75469 B 75469 ^ 

75459 P 75459 _ 

75449 T 75449 

75439 75369 



75359 Z 75279 v 

75269 Q 75079 _ _ _ 

75069 0060 

75059 N _ 

75049 M 

75039_?015_ 

75029 000 _ _ 

750 19 0J099 

7 5009 I 7.3.8 



14999 0000000000 _ 

14989 J 

74979 0030_ _ 

74969 0402 
74959 0400 

74949 0403 

74939 K 
74929 0007 

74919 L 

74909 0003 
74099 0009 

74889 0100 ______ _ _ 

74879 0010 
74869 00 T1 

74859 00 15 

14849 0404 

748 9 5043429400 _ _ 

74829 001 

74819 0405 _ _____ 

74809 0013 

74799 J5410000000__ 

74789 5110000000 _ _ 

74779 0020 ' __ ___ 

74769_001_4 

74759 0600 _ __ _ 

74749 0015 _ _ _ __ __ _ _ _ 

J4739_00J7 . 

74729 0016 
74719 0018 

74709 0019 __ ' _ 

74699 0041 

74689 SUMP 

74679 V 7.3.9 



74669 0051 
74659 0074 
7464-9 003 
74639 0050 
74629 E 
74619 003^ 
74609 ^8 1500000D 
74599 0200 
74589 0201 
74579 0210 
74569 0202 
74559 0310 
74549 0205 
74539 0212 

LOAD SUBROUTINES 

PDQ FIXED FMT SUBROUTNS 11/63 

PROCESSING COMPLETE 



7.3.10 
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© 



PROGRAM LISTING 

COMBUSTION PROCESS II 



PDQ FORTRAN C2 
START 



-6600 C SECOND PART OF COMBUSTION PROGRAM 

-6600 C NOVEMBER 6, 1964 

-6600 DIMENSION X(9,23) ,Y( 10,23) ,A( 1 ) ,C( 1) ,T ( 1 ) ,R(9, 1 0) ,P(20) ,V(9) 

-6600 READ 200, N,M 



-6636 C MULTIPLY X MATRIX BY Y MATRIX 

-6636 DO 22 1=1,9 

-6648 DO 22 K=1 ~10_ ^ 

-6660 R(I,K)«#.~ 



-6732 DO 22 J=1 ,23 

-6744 22 R(l,K)=R(i,K)+X(l,J)*Y(K,J) 



-7140 DO 23 1=1,20 



-7152 23 P(J)=X(8,I) 



-7284 DO 99 1=1,9 

-7296 DO 99 J=1 ,23 

-7308 99 X(l,J)=0.~ 



-7452 C GENERATE N MATRIX 



-7452 DO 3 1=1,9 



-7464 3 X(I,1)=R{I,1) 

-7572 DO 4 J=2,1'0 



-7584 4 X(1,J)=R(1 L J)_/X(1,1) 

-7752 DO 5 1=2,9* " _ 

-7764 5 X(I,2)=R(I,2)-X(1,2) *X(I,1 ) _ 

-7932 DO 55 J=3,10 " " 

-7944 55 X (2,J)=( R (2,J)-X(1,J)*X(2,1))/X(2,2) ^_ 

-8196 DO 44 l=2,9~ 7.3.II 



-8208 IK=I-1 

-8244 DO 33 J=l,9 

-8256 X(J,I)=R(J,I) 

-8400 DO 33 K=1 , IK 

-8412 33 X(J,I)=X(J,I)-X(K,I)*X(J,K) 

-8784 IF (X(|,|)) 34,56*34 " _ 

-8900 56 X(l,l)=i. " 

-8972 34 IM=I+1 

-9008 ^^-yjl'JL 

-9020 X(I,J)-R(I~J)/X(I,I) 

-9236 DO 35 K=1,IK 

-9248 35 X(l,J)=X(j,J)-X(K,J)*X(l,K)/X(l,l) 

-9692 44 CONTINUE 

-9728 C COMPUTE DELTA LOG T AND DELTA LOG P14-P20 

-9728 V(1)=X(9,10) . / . 

-9740 DO 6 J=2,9 _ ; 

-9752 K=10-J _ 

-9788_ L=J-1 _ . 

-9824 V(J)=X(K,10)_ 

-9896 DO 6_ 1-1 ~L _ 

-9908 IM=10-I 

-9944 6 V(J)=V{ J)- X(K,IM )*V(I) 

J0208 C COMPUTE DELTA* LOG P(l) AND S TORE IN R OW 9 OF X-MATRIX 

J0208 D O 10 J=1,20 

J0220 *[h{} m yi?» j wM- Y i}*»*L 1 

J0412_ DO J0 1=1,7 

^J0424 K= 1 0- 1 _~ 

J0460 10 X(9,J)=X(9,J)+Y(I,J)*V(K) _ _ 

J0760 C COMPUTE DELTA LOG P(N+1) STORE IN ROW 9 OF Y-MATRIX 

J0760 L=M 

J0772 M=M+13 

J0808 DO 1J_|=|,N 

J0820 11 Y(9,I)=L0G(P(I))*.434294+X(9,I) 



7.3,12 



J1024 DO 12 I=14,M 

J1036 12 Y(9,I)=L0G(P(I))*.434294+X(9,I) 

J1240 C COMPUTE DELTA LOG A(N+1 ) ,T(N+1 ) 

J 1240 A(1)=LOG(A(1))*.434294+V(2) 

J 1 288 T( 1 )=LOG (T ( 1 ) )*.~434294+V( 1 ) 

J1336 T(1)=EXP(T(1)/.434294) 

J1384 A(1)=EXP(A(1)/.~434294) 

J1432 PUNCH 102, T(1)" 

J1456 LL=101 

J1468 PUNCH 200,N,L,LL 

J1516 LL=201 

J1528 PUNCH 100,A(1),C(1),T(1),LL 

J1588 DO 7 1=1, N 

J 1600 LL=300+! 

J1636 7 PUNCH 101,Y(9,I),LL 

J1756 DO 8 1=14, M 

J 1768 LL=300+I 

J1804 8 PUNCH 101,Y(9,I),LL 

J1924 STOP 

J1936 200 FORMAT ( 13, I2,71X, 14) 

J1976 100 FORMAT (3E14.6,34X,|4) _ 

J2020 101 FORMAT (E14.6.62X, 14) 

J2054 102 FORMAT (F8.2)~ 

J2076 END 



T9999 SIN 
"19989 SINF 
T9979 COS 
79969 COSF 
79959 EXP 
79949 EXPF 
79939 LOG 
79929 LOGF 
79919 SQRT 



79909 


SQRTF 




79899 


ABS 




79889 


ABSF 




79879 


DRH 




79869 


DRHF 
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ATAN 
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79839 


X 


77779 


77769 
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75469 
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75459 


75449 


T 
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R 


74549 
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74349 
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V 


74259 


74249 


0200 
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74229 


M 
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0022 
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74199 
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000005 


90000 


74179 
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000 


— — - 


74159 


0023 




74149 


0099 




741 3 9 


0003 




74129 


0004 




74119 


0005 




74109 


0055 




74099 


001 




74089 


0044 




74079 


IK 





JL*2.3k 
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Tif069 0007_ ,; • f- '■ ■ 

74059 3033 _. ' : ->' 

0034 

'__ 74039 y*&>_ 

74029 5110000000 

74019 IM 

74009 0035 

73999 0006 _ 

739£9_0010__ ■_ _ _ 

73979 L . 

73969 0010 

73959 00 15 

73949 0011 

73939 5043429400 _ _ 

73929 0012 

73919 0102 

_ ][3909 LL 

73899 0107 __ 

73889 020 7 .V£ 

73879 0100 

_ 73869_0007 

73859_030_ 

73849 0101 

" 73839 0008 

18 

17 
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■ Fold under at dotted line . 
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• Fold under at dotted line . 
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7.5.1. 



• Fold under at dotted line. 



COEFFICIENTS OF THERMODYNAMIC DATA TABLES 



© 



TABLE II 



Coefficients and Constants for C p> Rj,, and S T Equations 



Tegp.°K 

Argon 

300-900 
900-3500 
1500-25C)0 
250O-35OO 

Br (gas) 

300-900 
900-1500 

1500-2500 
2500-3500 

Br 2 (gas) 

300-900 
900-1500 

i56o-256o 
2500-3500 

C (graphite) 

300-900 
900-1500 
1500-2500 
2500-3560 



C (gas) 



300-900 
900-1500 

1500-2500 
2500-3500 



U.96830 

1*.96830 
U.96830 
u.96830 



U.97787 
5.061*68 
3.9U566 



7.56238 
8.32286 
8.08510 



.1.37805 
1.7831*3 
3.69097 
5.26020 



5.011*1*6 
k.9hh01 
5.35597 
5.2U32U 



4x 30 3 



Y x 10 6 S x 30* 



0.00000 
0.00000 
0.00000 
0.00000 



- 0.01130 

- 0.67585 
1.624307 



0.00000 
0.00000 
0.00000 

0.00000 



- 0.16179 

1.0261*3 

- 0.58561 



No data presently available 



5.23251* 
1.U698U 
1.38077 



- 6.79762 

- 1.05952 

- 0.63636 



No data presently available 



3l*.3703£ 
5.76758 
2.06002 
0.1*5665 



0.17263 
0.081*99 
0.28602 
0.1*3913 



.30.56905 
. 3.01310 
. 0.55682 
. 0.01772 



0.22869 
0.09298 
0.12290 
0.201*30 



0.00000 
0.00000 
0.00000 
0.00000 



0.31111 
-0.30833 
0.06731 



3.05556 
0.27778 
0.101*90 



2.72222 
0.61111 
0.05206 
-0.00565 



-0.30278 
0.03333 
-0.0081*5 
-0.02211 



0.00000 
0.00000 
0.00000 
0.00000 



26.55610 
26.56008 
26.96793 



7.56818 
7.1*691*7 
7.66888 



92.03903 
91.18775 
90.1*3503 
89.30361* 



261.99066 
262.00556 
261.91501* 
261.881*75 



8.67690 
8.67698 
8.67698 
8.67701 



13.1*5302 
13.12831 
19.221*29 



3lu27127 
30.83361 
32.1*31*90 



5.381*93 
•30.92577 
.21.1*5576 
•31.09357 



9.23312 
9.57770 
8.38821* 
7.901*97 



ox 



7.6.1 







TABUS II (continued) 



Temp < °K 



CO (gas) 



30O-9OO 
900-1500 
1500-25QO 
2&0-3500 

CO2 (gas) 

300-900 
900-1500 
1500-2500 
2500-3500 

COS (gas) 

300-900 
900-1500 
1500-2500 
2500-3500 

CS>2 (gas) 

300-900 
900-1500 
1500-2500 
2500-3500 



H (gas) 



300-900 
900-1500 
1500-2500 
2500-3500 



Hg (gas) 



300-900 
900-1500 
1500-2500 
2500-3^00 



7.1*0128 
5.36912 
6.20126 
7.1*3752 



U.62716 
6.57128 
9.63218 
12.11555 



5.03988 
8.33898 
II.08096 
12.9U778 



6.31*585 
9.82266 
12.U3216 
13.723U8 



U.96810 
U.96810 
U.96810 
U.96810 



6.3UUU2 
7.6219U 
5.80531 
5.87503 



J x ID 3 



3.29913 
U.06773 
2.55577 
1.001*33 



18.19396 
U.0979U 
5.022U7 
1.886U1 



22.22610 
9.27936 
3o79860 
1.1*3330 



20.63U75 
7.8U001 
2.5180U 
0.86699 



0.00000 
0.00000 
0.00000 
0.00000 



3.0695U 
2.25177 
1.U70U3 
1.69951 



tx 10* 



7.16238 

- 1.803U5 

- 0.88801 

- 0.23500 



•1U.66U29 
. 5.82321 
. 1.767U0 
. 0.1*3936 



•22.1377U 
. 5.00821 
. 1.32033 
- 0.31585 



-2G.5396U 

- U.55571 

- 0.90U1U 

- 0.15699 



0.00000 
0.00000 
0.00000 
0.00000 



U.90655 
2.1*5762 

0.103U9 
0.32162 



Sx 10? 



•3.35833 
0.29722 
0.1121*9 
0.02038 



U.86389 
1.13611 
0.22673 
0.03831 



8.58889 
1.00556 
0.17053 
0.02768 



7.70000 
0.96111 
0.11890 
0.0175U 



0.00000 
0.00000 
0.00000 
0.00000 



2.76111 
-0.60833 

-0.01733 
0.02891* 



65.70215 
66.12651* 
65.78319 
65.0392U 



0.17012 

- 0.2l»259 

- 1.1*1059 

- 2.8951*6 



159.12375 
158.1*8017 
157.1*1*110 
156.32907 



317.51697 
316.79158 
315.8225U 
315.06061 



85.33693 
85.33691 
85.33690 
85.33690 



67.1*5133 
67.22029 
67.68U9U 
67.68U9U 



5.71*31*8 
15.67967 
11.03990 

3.68506 



19.89161* 
10.37873 
. 6.1*3266 

.21.15119 



20.89278 
5.0091*7 
•10.0311*9 
.21.11925 



15.36936 
. 1.60107 
.15.862U1 
-23.52007 



0.9138U 
0.91379 
0.91377 
0.91376 



- 5.6573U 
•11.72213 
. 1.80li05 
. 2.1*8129 



7.6.2 



TABID H (continued) 



Temp.°K 

l HBr (gas) 

300-900 
900-1500 
l500-2500 
2500-3500 



4' x lo- 



io 5 



H 2 (gas) 



300-900 
900-1500 

1500-2500 

2500-3500 



(gas) 



N (gas) 



300-900 

900-1500 
1500-2500 
2500-3500 



300-900 

900-1500 

i5oo-25bo 
25do-35bo 



h 2 "(«■*) 



300-900 
-900-1500 
l5oo-25oo 
25bo-35bo 



NO (gas) 



300-900 
900-1500 
1500-2500 
2500-3500 



7.1*031*7 
l*.82l*3l* 
5.61*971 



7.96561* 
7.16810 
1*.95233 
6.62109 



7.691*75 
5.19272 
6.06270 
9.11515 



1*.9681D 
1*.97213 
U.89391 
U.91189 



7.38656 
5.1*9191* 
5.95705 
7.20928 



7.83851* 
5.1*671*2 
6.61*779 
7.7459U 



- 2,78673 5.03202 -1.92778 

iu 73588 - 2.25298 0.1*3889 
2.93791 - 0.95U58 0.11886 

No data presently available 



1.18130 
1.7U67 
6.38682 
1U3561 



0.00779 
8.32851* 
7.0231*6 

3.21967 



0.00000 
0.01102 
0.13226 
0.11*696 



2.981*33 
3.1*1*012 
2.71090 
1.11*339 



5.0351*9 
i*.l*88l9 
2.210U* 
0.821*99 



5.26881 
1.531*76 
1.728H* 
0.961*69 



6.1*5881 
2.897H* 
2.30861 
0.71886 



0.00000 
0.00976 
0.07869 
0.09931* 



6.01631 
1.28357 
0.92791 
0.26995 



10.1*6357 
. 2.251*17 

- 0.77591 

- 0.19011* 



-2.22222 
-0.58611 
0.18121* 
0.08125 



-3.18056 
0.35556 
0.28578 
0.06309 



0.00000 
-0.00278 
0.01568 
0.0201*1* 



-2.61111 
0.16667 
0.11599 
0.0231*5 



•5.20833 
0.1*2222 
0.09955 
0.01657 



25.26861 
25.93819 
25.61*681* 



11.33351 
11.56568 

12.30309 
11.22879 



162.2759U 
162.81*1*91 
162.1*31*140 
160.581*88 



113.1*2918 
U3.1*28U 
113.1*601*0 
113.1*3056 



1.651*02 
2.07553 
1.86228 
1.10708 



23.31*538 
23.78872 
23.32637 
22.66911 



5.90995 
19.08726 
U*.6l765 



. 0.13970 
3.76566 

15.80923 
5.76571* 



5.061*17 
17.51098 
12.52268 
. 5.6581*6 



8.30830 
8.28757 
8.7231*3 
8.58571 



1*.33152 
13.71896 
11.06826 

3.6151*3 



6.7661*3 
18.10630 
11.59108 

5.06361 



TABLE II (continued) 



Temp.°K, 



N0 2 (gas) 



(gas) 



300-900 
900-3500 

1500-2500 

2500-3500 



300-900 
900-1500 
1500-2500 
25GO-3500 



2 (gas) 



300-900 
900-1500 
1500-2500 
2500-3500 



CH (gas) 



S (gas) 



300-900 
900-1500 
1500-2500 
2500-3^00 



300-900 
-900-1500 
l50O-2$OO 
2500-3500 



S 2 (gas) 



300-900 

90o-i5oo 
1500-2500 

2500-3500 



6,07017 
6.U3739 
10.10326 
12.06Ul5 



5.821A7 
5.1825U 
5.01322 
5.15709 



6.86951; 
5.6U527 
7.61206 
7.82021 



7.58972 
7.20029 
5.38982 
5.89037 



6.01259 
5.90721 
5.71818 
5.20717 



5.70207 
7.77137 
8.U8089 
8.72567 



6 x 10- 



10.18U7U 
11.25150 
3.8U337 
1.3571*5 



2.93332 
0.351^3 
0.01167 
0.17073 



0.97686 
U.87180 
0.9W&7 
0.58387 



2.19358 
1.33i£7 
2.U3528 
1.88039 



1.111*69 
1.39923 
0.9U930 
0.37068 



9.85685 
2.0032U 
0.53176 
0.211*U3 



Yjl 10 



- 2.27786 

- 6.39881 

- 1.35869 

- 0.302U7 



3.68702 
0.21798 
0.0!l«8U 
0.0U35U 



6.00619 
2.79798 
0.16012 
0.031U6 



2.50U6U 
2.02286 

o.6ioia 

0.140655 



0.39298 
0.770U8 
0.U1903 
0.200^2 



-11.21*595 

- 1.20259 

- 0.17635 

- 0.03888 



S x 10' 



-l.li36ll 
1.33356 
0.1769U 
0.02667 



-1.60000 
-0.0U722 

0.00591 
-0.00120 



-3.63056 
0.61667 
0.02108 

-0.01127 



-0.56667 
-0.55833 
0.05829 
0.03351 



0.65833 
-O.U1I67 
-0.05130 
-0.0237U 



U.569141* 
0.26389 
0.023U5 
0.00355 



11.76362 

11.5U805 
10.173U9 

9.00610 



59.97630 
60.09763 
60.16089 
60.06363 



2.05967 
2.23511 
1.1*8991 
1.141325 



1*3.90179 
Uk.01996 
UU.67687 
lOt.33976 



l6U.700l|0 
16U.75922 
16U.81760 
165.15635 



228.792li8 
228.37683 
228.11776 
227.97571 



19.83972 
17.377U6 
. 2.69188 
-llt.33726 



6.00725 
9.08008 
10.00671 
9.13323 



9.92168 
15.51879 
U.7293U 
3.57106 



1.19689 
3.26569 
13.12031* 
10.08323 



6.171*1*7 
6.87100 
7.87223 
10.9631*0 



19.5U3U9 
9.5ll*70 
5.61*902 

U. 20115 



TABI3 II (continued) 



(gas) 



ID 3 X * ID 6 $ x ID? 



(gSS) 300-900 5.9530U Hw686l2 -9.7131*3 2.01*722 30.81721 J^^Sg 

W-1SX> 7?2609 9.I.22U2 -5.50357 1.1*91* 30.29^0 11.23559 

iSn £fon 11 10«>!i6 2.92319 - 1.03367 0.13607 29.11630 - 6.13302 

2§MSE SlO 1.01W1 -O^&S 0.019U1 28.23600 -lU.97395 



300-900 U.53177 33.1531*2 -28<,72571 ' 9.21910* 9*10707 26.83W5 

9m-l5b0 9!87962 15.56309 - 9.25952 1.97222 7.87669 0.16531 

lISS 15 23516 U 60266 - 1.71281* 0.22566 5.89l*5U "29.08528 

25^-3500 17 71931 l.**32l*5 - 0.35721 0.03160 Wi25« -1*3.8211*9 



© 







Tewp«°S 



TABUS HI 

Coefficients and Constanta for -Aq/foX Equations 

a b x 10^ cxlO 6 d x 10? 



Br 2 (gas) 



300-900 
900-1500 
1500-2500 
^00-3500 



00 (gas) 



300*900 
900-3500 
1500-2^00 
2500-3500 



00 2 (gee) 



300-900 
900-1500 

i5oo-25oo 

2500-3500 



COS (gas) 



300-900 
900-1500 

l5oo-25oo 
25oo-35oo 



CSj (gas) 



300-900 
900-1500 

1500-2500 
2500-3500 



1.201O6 
0.90905 
-0*09751 



1.7298k 
2.39396 
1.99668 
1.1*9090 



6.05670 
k.39692 
2.79290 
1.73197 



5.91*368 
3.87205 
2.1*1861 
1.33831* 



5.38117 
3.1*9011 
2.093U2 
0.97325 



- 1.32211 -1.08593 -0.30612 

- 0.70996 0.52206 -0.11252 
0.1*791*0 -0.08971 0.00371* 

lo data presently available 



0.01*861 
1.09122 
0.71791* 
0j*06l3 



6.09715 
2.9l*9ll* 
1.31*150 
0.67102 



6.65107 
2.751*50 
1.2691*8 
0.60733 



5.79611 
2.65528 
1.18320 
0.51515 



-0.5UU58 
0.323U7 
0.16708 
0.08099 



3.731*96 
1.031*29 
0.31209 
0.12157 



k.30l*l6 
0.99026 
0.30988 
0.12817 



3.35175 
1.00703 
0.3128U 
0.13U55 



0.20827 
•0.03911* 
-O .011*1*7 
-0.00550 



-1.02738 
-O.15061 
-0.02810 
-0.00790 



-1.21189 
-0.11*607 
-0.02823 
-0.0091*0 



-0.81596 
-0.15236 
-0.02893 
-0.01096 



22.91800 
22.97167 
23.281B0 



128.9538* 
128.80881* 
128.96790 
129.27807 



192.11031 
192.1*1*760 
193.05356 
193.68752 



16U.82172 
165.21*372 
165.78225 
166.1*1*810 



137.81513 
138.21*681* 
138.71*767 
139.U5676 



7,6.6 



TABUS III (continued) 



% (gas) 



Tnp.°K 



300-900 
900-1560 
1500-2^00 
2500-3500 



BBr (gas) 



300-900 
900-1500 

i5oo-25oo 

2500-3500 



HaO (c*») 



300-900 
900-1500 

1500-2500 
2500-3500 



%•■ (|M) 



■2 (f») 



WO (gaa) 



100-900 
900-1500 
1500-2500 
250O-350O 



300-900 
900-1500 
1500-2500 
2560-3500 



300-900 
90O-150O 

i5oo-25bo 
25bo-35oo 



1.80740 
i.iOi55 
2.07869 
2.04360 



1.27940 
2.62092 



3*92280 

luooodo 
5.O3059 
4.26326 



Iul53b8 
5.359W 
4.82659 
3.03344 



1*28299 

2.2401*3 
1.92766 
1*31563 



1.4861*9 
2.358« 
1.61t012 
1.16897 



b x 10- 



-0.77231 
0.56655 
-0.36996 
- O.I|276o 



0.69836 

- 1.35657 

- 0.32578 
■0 data prat 



• 0.4I1081 

- 0.52733 

- 1*60988 

- 1.15897 



- 0.27850 

- 2*1*753 

- 2.00597 

- 0.90334 



O.75087 
0.87109 
0*61552 
0.21378 



0.52891 
1.22114 
0.52574 
0.21358 



c x 10 



0.82300 
-0*10223 
0.01736 
0.05395 



d x 10 ; 



-0.3U735 
0.07553 
0.00218 

-O.0036U 



-0.87118 
0.55007 
0.06189 
intly avallabl* 



0.28166 
-0.091*00 
-0.0062(9 



-0.26532 
-0.22087 
0.28738 
0.16912 



-1.1*929 

0*61519 
0.k5752 
0.15U20 



-1*00915 
0*21858 

0.12924 
0.01196 



-1.13667 
0.41630 
0.11446 
0.02253 



0*07828 

0.06779 
-0.02206 

-0*01037 



0*U8294 
-0.06255 
-0.04241 
-0.01092 



0.32848 
-0.02167 
-0.01065 

0.00219 



0.45393 
-0.059lil 
-0.00981 

0.00033 



51.91*214 
52*05838 
51*72101 
51*821*55 



43*58968 
43*25493 
43.60677 



110*36133 
110.33572 

109*96629 
110.45794 



87.10401 
86.84728 
87*08322 
88.10437 



113.32404 
113.11086 

113.25067 
113*60060 



75.5U06 
75.34848 
75.62922 
75*89597 



7.6.7 



6 



© 



a 



TABLE III (continued) 



Temp.°K 



02 (g*») 



CH (gas) 



&2 (■*•) 



300-900 

9oo-i5bo 
i5oo-25oo 
256o-356o 



300-900 
900-1^00 

1500-2500 
25oo-35oo 



300-900 

900-1500 
i5oo-25oo 

2500-3500 



300*900 

9oo-i5oo 
i5oo-^oo 
25oo-35oo 



SC^ (gas) 



SO3 (gas) 



300-900 
900-1500 

1500-2500 
2500-3500 



300-900 
900-1500 
l5oo-25oo 
2560-3500 



5.03721* 
1*.1*781*1* 
2.1*2399 
1.59110 



2.1t0502 
2.37563 
1.21U93 
1.251*98 



1.61169 
1.1*81*61* 
2.3101*7 
2.13098 



3.18182 

2.031*1*9 
1.1*8721 
0.81*975 



5.89177 
U.19986 
2.331*1*1* 
1.1*7201 



9.53787 
5.821*71 
2.77901* 
1.1*8906 



b x 10- 



U.03857 
3.01193 
1.87920 
0*78095 



1.23028 
l.ljOl*01 
0.21*31*0 
0,23281 



0.18612 
0.21*658 
0.61566 
0.51607 



3.01*093 
1.20812 
0.6Uli9 
0.21*01*8 



5.1*5159 
2.90101 
0.98020 
0.1*3558 



•10.83605 
■ fc.535H* 
. I.li0570 
- 0.58201 



cxlD' 



1.61896 
1.1U807 
0.62917 
O.ll*6ol* 



0.2291*1* 
0.51*21*1* 
0.02188 
0.00933 



0.19833 
-0.30271* 
0.09990 
0.07550 



1.751*51 
O.U6019 
0.17015 

0.07376 



2.79992 

1.12550 
0.23869 
0.08U95 



6.60773 
1.79207 
0.35612 
0.115U5 



d x 10 : 



-O.22150 
-0.17962 
-0.07520 
-0.00l*3l* 



0.05U16 

-0.0891*6 

-0.00117 
0.00111 



-0.13000 
0.06U5D 
-0.00659 
-0.001*37 



-0.1*0920 
-0.06881* 
-0.01596 
-0.0061*2 



-0.57729 
-0.17682 
-0.02209 
-0.00573 



-1.68085 
-0.28375 
-0.03261 
-0.007U2 



111.5191*6 
111. 7U950 
112.52110 
112.99561 



59.321*36 
59.35819 
59.79681* 
59.73753 



51.03079 
51.03237 
50.73363 
50.851*33 



50.62665 
56.89561 
51.081*13 
51.1*9653 



127.73151* 
128.11*633 
128.83220 
129.31*775 



168.83660 
169.66850 
170.79080 
171.55359 



7.6.8 



TABLE IV 



Br 2 (gas) 



00 (gas) 



Temp.^E 



300-900 
900-1500 

l5bo-25oo 

25bO-3500 



300-900 
900-1500 

i5oo-25bo 

2500-3500 



Coefficients and Constants for Log K Equations 
A B x ID 3 P x ID 6 A x 10? 



C0 2 (gas) 



300-900 
900-1500 

i5bo-25bo 

2500-3500 



COS (gas) 



300-900 
900-2500 

i5oo-25oo 
25oo-35bo 



CS2 (gas) 



300-900 
900-1500 

i5bo-25bo 
25oo-35oo 



-1.20U36 
-0.90905 
0.09751 



.1.72981* 
•2.39396 
•1.99668 
.1.1*9090 



-6.05670 
-t.39692 
-2.79290 
-1.73197 



-5.91*368 
-3.87205 
-2.1*1861 

-1.33831* 



-5.38117 
-3.1*9011 
-2.0931*2 
-0.97325 



0.57U8 -0.23581 0.01*1*31 
0.30833 -0.11336 0.01629 
-0.20820 0.0191*8 -0.00051* 
Ho data presently available 



M 



-0.02111 
0.1*7391 
0.31180 
0.17638 



2.61*796 

1.28079 
0.58261 
0.2911*2 



2.88983 
1.19627 
0.55133 
0.26376 



2.51722 
1.15317 
0.51386 

0.22373 



0.11825 
-0.07021* 
-0.03628 
-0.01759 



-0.81101* 
-0.221*59 
-0.06777 
-0.0261*0 



-0.931*61* 
-0.21563 
-0.06729 
-0.02783 



-0.72782 
-0.21867 
-0.06793 
-0.02922 



-0.03015 
0.00567 
0.00210 
0.00080 



0.11*873 
0.02180 
0.001*07 
0.00 111* 



0.1751*1* 
0.02115 
0.001*09 
0.00136 



0.11812 
0.02206 
0.001*19 
0.00159 



2.23727 
2.97500 
5.72736 



1.32l*lJB 
0.38885 
0.71*022 
2.27070 



2.33l<06 
1.881*10 
6.39973 
9.5981*1 



2.1*6722 
. 2.80267 
. 6.88230 
-10.15361 



0.97930 
. 3.9301*3 
• 7.83135 
.11.23680 



9.95316 
9.9761*7 
10.11116 



56.0039U 
55.9k)97 
56.01005 
56.110*75 



83.1*321*5 
83.57893 
83.81*207 
81*.U7l*2 



71.58117 
71.761*1*1* 
71.99832 
72.2871*9 



59.85235 
60.03981* 
60.25735 
60.56530 



7.6.9 



C 



o 



H2 (gas) 



300-900 
900-1500 
3500-2500 
2500-3500 



HBr (gas) 



300-500 
900-1500 
1500-2500 
2^00-3500 



EjO (gas) 



300-900 
900-X500 

i5oo-25oo 
2500-3500 



H 2 S 



(gas) 



No (gas) 



HO (gas) 



300-900 
900-1500 
1500-2500 
2500-3500 



300-900 
900-1500 
1500-2500 
2500-3500 



300-900 
900-1500 

1500-2500 
2500-3500 



•1.8071*0 
.1.161*55 
.2.07869 
.2.014360 



.1.279140 
-2.62092 
•1.61*21*9 



-3.92280 
-1*.00080 
-5.03059 
-U. 26326 



-l*.l53l*8 
-5.3591*8 
-U.82659 
-3.033W* 



•1.28299 
-2.21*010* 
.1,92766 
•1.31563 



.l.lf861*9 
.2.35865 
.1.61*012 

•1.16897 



TABLE IV (continued) 
B x lp3 r x IP 6 a x 10? 



0.3351a 
-0.21*605 
0.16067 

0.18570 



-0.30329 
0.58915 
0.3ial*9 



-0.17871 
0.08951 
-0.00377 
-0.01171 



0.18918 
-0.1191*5 
-0.0131*1* 



0.05628 
-0.01108 
-0.00032 

0.00053 



-0.01*077 
0.01361 
0.00091* 



No data presently available 



0.1911*1* 
0.22902 
0.69916 
0.56333 



0.12095 
1.06295 
0.87118 
0.39232 



-0.32610 
0.37831 
0.26732 
0.09281* 



-0.22970 
0.53033 
0.22832 

0.09275 



0.05761 
0.01*796 
-0.0621*0 
-0.03672 



0.21*956 
-0.13359 
-0.09935 
-0.0331*8 



0.21913 
-0.01*71*6 
-0.02806 
-0.00260 



0.21*682 
-0.0901*0 
-0.021*85 
-0.001*89 



-0.01133 
-0.00981 
0.00319 
0.00150 



-0.06991 
0.00906 
0.00611* 
0.00158 



-0.01*755 
0.00311* 

0.00151* 
-0.00032 



-0.06571 

0.00860 

0.0011*2 

-0.00005 



- 0.05198 
.1.65659 
0.90789 
0.71*1*65 



O.89267 
2.61*026 
0.09371 



0.75972 

0.97539 
3.85222 

1.511*97 



1.96065 
5.05225 
3.51185 
1.91571* 



2.12757 
0.31*882 
0.55680 
2.39115 



1.001*21 
1.18577 
0.81*787 
2.25801 



M 



22.55819 
22.60867 
22.1*6215 
22.50712 



18.93081* 
18.78538 
18.93818 



U7.92932 
1*7.91820 
1*7.75775 
1*7.97128 



37.82879 
37.71730 
37.81976 
38.29799 



1*9.21601 
1*9.1231*2 
l*9.l£l*U* 
1*9.33612 



32.791*01* 
32.7231*3 
32.81*535 
32.96120 



Terap.°K 

NO2 (gas) 

300-900 
900-1560 
1560-2500 
2500-3500 

O2 (gas) 

300-900 

900-1500 
1560-2500 
2500-3500 

OH (gas) 

300-900 
900-1500 
1500-2500 
2500-3560 

S 2 (gas) 

300-900 
900-1560 
l5oo-25oo 
2500-3500 

502 (gas) 

300-900 
900-1560 
1500-2500 
2500-3500 

503 (gas) 

300-900 
900-1560 
1560-2500 
2500-3500 



(continued) 

p x 10 6 & x 10 9 L M 



A 



-5.30721* 
-l*.li 781*1* 
-2.1.2399 
-1.59110 



-2.1*0562 
-2.37503 
-1.231*93 
-1.251*98 



-1.61169 
-1.1*81*61* 
-2.33X>1*7 
-2.33098 



-3.18182 
-2.031*1*9 
-1.148721 
-0.81*975 



-5.69177 
-li.39986 

-2.331*1*1* 
-1.1*7201 



-9.53787 
-5.821*71 
-2.77901* 
-1.1(8906 



TABLE IV 
B x 10 3 



1.75393 
1.30807 
0.81613 
0.33916 



0.53U31 
0.60975 
0.10571 
0.10111 



0.08083 
-0,30709 
0.26738 
0.22U12 



1.32066 
0.521*68 
0.26557 
0.101*1*1* 



2.36760 
1.25989 
0.1*2570 
0.3B917 



1.70603 
1.96959 
0.610U9 
0.25276 



-0.35155 
-0.21*930 
-0.13662 
-0.03171 



-0 .01*982 - 
-0.11779 
-0.001*75 
-0.00203 



-0.01*307 
0.06571* 
-0.02169 
-0.01639 



-0.38099 
-0.09993 
-0.03695 
-0.01602 



-0.60799 
-0.2l*i*iiO 
-0.05183 
-0.0181*5 



-1.1*31*85 
-0.3893J* 
-0.07603 
-0.02507 



0.03212 
0.02600 
0.01089 
0.00063 



-0.00781* 
0.01295 
0.00017 

-0.00017 



0.01882 
-0.00931 
0.00095 
0.00063 



0.05921* 
0.00997 
0.00230 
0.00093 



0.08357 
0.02560 
0.00320 
0.00083 



0.21*333 
0.01*108 
0.001*72 

0.00307 



2.1993U 

- 0.03725 

- 2.1*9391 

- 6.31i*51* 



0.58712 
0.1*51*22 

- 2.81251* 

- 2.6661*9 



- 0.15157 

- o.li262l 
1.88357 
1.33278 



2.951*11* 

- 0.01*026 

- 1.56036 

- 3.501*71 



3.26278 

- 1.19090 

- 6.1*2058 

- 9.02099 



U.71887 
- lt.88888 
-13.1*3033 

-17.31388 



1*8.1*3229 
1*8.53219 
1*8.86729 
1*9.07337 



25.761*21* 
25.77891* 
25.9691*1* 
25.91*368 



22.16239 
22.16308 

22.03331* 
22.08576 



21.98688 
22.1031*2 
22.18556 
22.361*66 



55.1*7311 
55.65325 
55.95112 
56.17502 



73.321*81 
73.68610 
7U.17356 

7l*.5ol*78 



